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INTRODUCTION 


Aircraft propulsion currently is almost completely 
dominated by gas turbines with essentially isobaric combus- 
tion so that the thermodynamic cycle approximates the 
classical Brayton cycle with isentropic compression and 
expansion, and isobaric heat exchange. 

It has been known virtually since the beginning of gas 
turbine development that significant gains in efficiency 
could be achieved if energy release could be accomplished 
at constant, volume, rather than at constant pressure. The 
constant volume or Humphrey cycle is compared with the 
constant pressure Brayton cycle in Figure 1 in pressure- 
volume and temperature-entropy thermodynamic diagrams 



Figure 1. Comparison of the Brayton and Humphrey Cycles 


The ideal gas thermodynamic efficiency of the closed 
Humphrey cycle is given by: 

(T 3 /T 2 ) l/y - 1 

= 1 - ( 1 ) 

T3/T2 - 1 
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and that of the Brayton cycle by 

= 1 - 1/R^ (2) 

Here p, v, and T, following convention, denote the pres- 
sure, volume, and temperature respectively and 

R - pressure ratio = P 2 /P 1 

7 - isentropic exponent 

r - 1 - 1/7 


Numerical subscripts denote the thermodynamic states in the 
cycle. Comparisons of the performance of the two cycles 
are shown in Figure 2. The displayed results indicate the 
significant superiority of the Humphrey cycle to the 
Brayton cycle even at modest temperature ratios. 

Numerous attempts have been made to achieve constant 
volume combustion, but the complications of valves, such as 
those used by Holzwarth, did not result in a practical 
engine. Recent improvements in gas turbine performance 
have diminished somewhat the interest in the development of 
alternative engine cycles. However, improvements in 
performance have been gained at the expense of the cost and 
complications of multistage axial compressors. Such 
compressors are economically unacceptable for smaller gas 
turbine engines which might find applications in light 
aircraft or cruise missile. 

Presented here is the summary of the second phase of a 
study aimed at demonstrating the feasibility of employing 
transient transverse detonation waves to augment the 
performance of gas turbines using low pressure ratio com- 
pressors. The ultimate objective is to replace multistage 
axial compressors in smaller gas turbines with single stage 
centrifugal units augmented by an array of detonation 
ducts. Alternately, specific fuel consumption and thrust 
coefficients of low-to-medium pressure ratio gas turbines 
could be enhanced through the use of detonation ducts. In 
the initial study (Reference 1) performance of pure jet 
engines with low pressure ratio compressors augmented by 
detonation ducts was studied parametrically , and it was 
shown that significant reductions in specific fuel con- 
sumption could be achieved concurrently with greatly 
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enhanced thrust coefficients. The main objective of the 
present study is to analyze the gas dynamics of the detona- 
tion duct-combustor interactions with the aim of determin- 
ing the realistic bounds of performance in actual engines. 

In the interest of clarity of presentation of the main 
objective of the program, details of the detonation duct 
gas dynamics have been placed in Appendix A, the supporting 
experimental work is in Appendix B, and the listings of 
computer codes developed specifically for this study are in 
Appendix C. 
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Figure 2. Comparison of Humphrey and Brayton 
Gas Turbine Cycles 




THEORETICAL CONSIDERATIONS 


The theoretical upper limit on performance could be 
achieved if the whole compressor output was processed by 
detonation waves. Lower limit corresponds to 100% of the 
compressor flow going through the combustor, which is the 
case with conventional gas turbines. It is of some 
interest to establish the upper limit as a measure of the 
potential performance improvements. 


Detonation Wave Engines 

In this case, the heat addition is by means of detona- 
tion waves. It is known (e.g., Stanyukovich^ ) that a 
Chapman-Jouguet detonation is essentially equivalent to a 
constant volume heat release. In the strong wave limit, 
the speed of the detonation wave, D, is given by: 


d 2 = 2(y2 - 1)Q 


(3) 

with Q being the energy release 
detonated mixture. The pressure 

relations are: 

per unit mass 
temperature and 

of the 
density 

PD = pd 2/ ( y + 1) 


(4a) 

Td = 2 y (Q/Cv)/( y + 1) 


(4b) 

pd = p(y + i)/y 


(4c) 


with the subscript D denoting conditions downstream of the 
detonation wave. Because of the assumption of the Chapman- 
Jouguet state, the gas behind the detonation wave is at 
Mach 1. Assumption of isentropic compression to stagnation 
conditions and ideal compressor and turbine performance 
results in the relation for engine efficiency which is 
given by: 

e = yD2(l - Rd-^D) - y/H(R^ - 1) (5) 
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with 


H = Q/Cv Ti 


Maximum efficiency is then achieved at: 


Ropt 


y / y - 1 

^ ^ (2H _D 

yy-i\ y-1 

D 


1 / 


y 

yf+ f 

D 


(6) 


with >'d being the isentropic process exponent of the 
detonation products. Variations of efficiency with pres- 
sure ratios for a range of values of H is presented in 
Figure 3. For hydrocarbon fuels, H ranges from 5 to 10. 
It should be noted that the efficiency of about 53%, which 
is reached by the proposed cycle at a compressor pressure 
ratio of 5, is matched by a standard Brayton cycle at a 
pressure ratio of about 14. Alternatively, augmentation of 
a R = 5 Brayton cycle with detonation wave heat release 
increases the efficiency from 37% to 53%. 
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Figure 3. Performance of Augmented and 
Brayton Gas Turbines 
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A theoretical detonation wave engine considered above 
cannot be realized because of practical considerations such 
as complete detonation of the whole flow, starting prob- 
lems, unacceptably high temperatures and maintenance of 
flow through the engine. An engine which could be made to 
work is shown in Figure 4 in pure jet engine form. 



DETONATION ® EXPANSION 
CHAMBERS 

(N) DENOTE STATION NUMBERS 

Figure 4 . Schematic Representation of a Detonation 
Wave Augmented Jet Engine 


The station numbers are given for easy reference. The 
detonation duct could be a rectangular channel such as is 
shown in Figure 5. An array of injectors sends pulses of 
gaseous, or liquid, fuel into the duct where the mixture of 
air and fuel is detonated by an array of igniters fired in 
sequence with the injectors, or operated continuously. 
Maximum rate of energy release is obtained when the flow in 
the channel is at Mach 1. The rapid area increase at the 
upstream channel is necessary for the attenuation of shock 
waves. Area increase at the downstream end may not be 
necessary with ejector mixing-wave compression mode of 
operation. Preliminary experimental studies using a 
detonation duct with spark plug igniters are described in 
Appendix B. 
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Figure 5. Schematic Representation of 
the Detonation Duct 


The range of operating conditions for the engine is 
bounded by no detonation duct flow and the maximum fre- 
quency cycling of the ducts. In the previous study, i1: was 
assumed that the combustor and detonation duct flows were 
mixed in an ejector. Since the cycling frequency of the 
duct was not determined at that time, the fraction of flow 
through the combustor was used as a free parameter. The 
studies were confined to a pure jet engine and representa- 
tive performance parameters for a pressure ratio of 4 are 
shown in Figure 6. The limit of =1 represents no 
augmentation by detonation ducts . Even though the concept 
is aimed primarily at the use of very low pressure ratio 
compressors, it is interesting to note in Figure 7 that 
significant gains are possible with moderate pressure 
ratios . 

In the present study, the pure ejector mixing model 
was rejected in favor of a mixed ejector-wave compression 
mixing of combustor and detonation duct flows. The more 
efficient mixing offset the lower than expected detonation 
wave compression which was measured in a preliminary 
experimental program. 
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Figure 7. Influence of Gas Dynamic Models 
on Engine Performance 
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BASIC CALCULATIONS 


All of the physical quantities in the detonation wave 
augmentated gas turbines are scaled on the detonation duct 
parameters evaluated prior to the detonations and denoted 
by the subscript o. The natural scaling quantities are the 
detonation volume acoustic velocity, ao, detonation volume 
cross-sectional area, and length, Aq and Lq detonation 
volume mass PoAqLo/ and internal energy CvTq. In the 
following analyses, these quantities appear: 

Time 


d — tao/Lo 


Area and length 


A 


Mass 


m 


Energy 


E 


A/Aq and x/Lo 


Po^oAq 


E/ Poi'oAocv'i’o 


All thermodynamic quantities and velocities scale on 
corresponding detonation volume parameters. 


Baseline Gas Turbine 


Performance of detonation wave augmented gas turbines 
is compared with a corresponding conventional turbine with 
pressure ratio R, compressor and turbine efficiencies ec 
and eT respectively, and a maximum to ambient temperature 
ratio <i> = Tmax/Ta- The efficiency of such a turbine is 
given by: 

©c eT ^ “ 1 

r) = ^ (7) 

e„( (b - 1 w rpi _ 1 ^ 1 
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Compressor 


The total temperature of a compressor with a pressure 
ratio R and an isentropic efficiency e^ is: 


TT/Ta = 1 + (r’" - l)/ec ( 8 ) 


The ambient temperature Ta is related to Tq through the 
total temperature and isentropic relations. Compressor 
work, nondimensionalized by the detonated volume energy, is 
given by: 


y(md + mb) [1 + ( y - l)Mo2/2] 
1 + ec/R^ - 1) 


(9) 


Detonation Duct 


The dimensionless heat released in the detonation 
volume is: 


Q = (Pd/^o)(TD/Td - 1) = (Pd/^o) [ (CD/Cd)2 - 1] (10) 


with the subscripts d, D denoting the conditions prior to 
and subsequent to detonation. In the calculations, the 
acoustic velocity ratio is taken as 2.5 and therefore the 
energy release in the three sequenced firings is: 


Ql = 5.25 Q 2 = 6.40 Q 3 = 8.24 

Detailed calculations of the gas dynamics of the duct and 
its output for various modes of operation are given in 
Appendix A. 


Combustor 


Heat is assumed to be added isobarically during the 
period A 6 . The dimensionless mass is: 


m — Mb(Pb/^o) (*-b/Co) (Ab/Ao)A(9 


( 11 ) 
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Initial combustor and detonation duct quantities are 
related through isentropic flow assumptions and Mach 

numbers and Mq. The temperature at the completion of 
combustion is: 

Tbf/Tb = 1 + Qc/Cp Tb = 1 + Q* (12) 

with Q* being the heat release parameter. The total 

dimensionless heat added is: 

Qb = yQ*(Tb/To)'^ (Ab/Ao)Mb Aa (13) 

with 

A = (3 y - i)/2 ( y - 1) 

Taking 


ec = 0.85 Mb = 0.25 Q* = 2 R = 5 
y = 1.32 

the maximum combustor temperature is calculated to be 
Tbf/Ta = 4.82 which corresponds to a turbine inlet tempera- 
ture of about 2000°F. In order to maintain a constant Mach 
number in the combustor under the condition of the above 
parameters, the combustor flow area must increase iDy a 
factor of 3 1/2 between the entry and exit. 


Mixing and Wave Compression 

The output of the combustor is steady while that of 
the detonation duct alternates volumes of fairly steady 
flow with sharp pulses of high pressure and temperature 
gas. When the detonation duct output is almost steady, the 
streams are mixed as if in an ejector with moderate pres- 
sure changes. At exit from the detonation duct, the 
pressure pulses expand suddenly and because of their much 
higher velocity and pressure, immediately compress the 
preceding volumes of mixed gases. In the initial phase of 
this process, the overall volume remains virtually 
unchanged while the flow equilibrates to an essentially 
uniform pressure and velocity. In the proposed system, the 
flows of 2 detonation ducts are mixed with the outflow from 
1 com.bustor. The arrangement is something like that shown 
in Figure 8. Further equilibration of the flow could be 
achieved in the optional spiral collector. Without the 
spiral collector, the turbine would be operated in purely 
blov/down mode. 
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TURBINE NOZZLES 


ALTERNATIVE BURNER AND 
DETONATION DUCT PORTS 



Figure 8 . Collector and Gas Dynamic Augmentor 



The mixing-wave compression section mixes the com- 
bustor and detonation duct flows and compresses them by the 
detonation wave moving downstream. The mixing is accomp- 
lished by momentum transfer as in ejectors. The high 
pressure detonation wave expands into the collector section 
and transfers energy to the mixed gases by doing expansion 
work on them. Curvature of the collector is considered 
necessary because in a real engine, this would add some 
length to ensure good mixing and would increase pressure by 
reflections of the detonation waves. A schematic represen- 
tation of the collector section is shown in Figure 9. 



Figure 9. Mixing-Wave Compression Section 
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The conditions at the end of mixing of the detonation 
duct and combustor streams are determined from the solution 
of the continuity, momentum, and energy equations. With 
all quantities normalized by the initial detonation duct 
values, the governing equations are: 

Continuity 


(Piti/Tm) A-m - S^i ~ ^m 


(14) 


Momentum 


E(Pi AiA0 + Ui mi) = Pm Am^^+ Um Mm (15) 

Energy 

E^i[Ti + (y - l)Ui2/2] = [Tm + (Y - l)Um2/2]Mm 

(16) 

with the definitions 

Sm (Tm/Um + Um)Mm 


and 


Se = [Tm + ( y - l)Um2/2]Mm 

Subscripts i denote the individual volumes of gases enter- 
ing the mixing section. The solution for the mixed velo- 
city is: 


Sm 

' / Sm \ 2 

2Se 

(3 - y )Mm 

\ (3 - y )Mm / 

(3 - y )Mm 


(17) 

The remaining quantities are then evaluated by simple 
substitution. 

When the high pressure volume of the shock wave com- 
pressed flow exits the detonation duct, it undergoes a 
sudden expansion into an area Am- Accompanying this sudden 
expansion is an entrepy increase or the so called "shock 
loss" which is accounted for here using the relations given 
by Rudinger. Immediately after expansion, the gas 
exchanges energy with the mixed gases by decelerating and 
expanding 
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to some common pressure Pe in a flow area Ae« In the first 
approximation, the equilibration process may be considered 
to be isentropic. 

The continuity equation is written as: 




Ai 



1/y 

A^i 


(18) 


The statement of constant total energy is: 


2 ; Ti 

y(y - ( 


/Pe \ y-l/K 


U±2 i = Mm Ue2 


(19) 


When the wave compression exit area Ae is specified, the 2 
equations are solved by iteration for Pe and Ue* This 
determines the thermodynamic state points of the wave 
compressed gas. The parameters of the mixing wave compres- 
sion process are: 


Ab/Ad “ combustor flow area 
Ae/Ad “ wave compression exit area 


Am/Ad “ mixing region area 


Other parameters which bear directly on the problem are the 
combustor conditions and single or multiple firings of the 
detonation duct. A representative example of the output of 
the mixing-compression section is shown in Figure 10. 
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Figure 10. Gas States at Collector Exit 
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Parametric Studies 


Extensive parametric studies were performed to deter- 
mine the efficiency of detonation wave enhancement of the 
performance of low pressure ratio gas turbines. The basic 
engine parameters were: 


R 

= pressure ratio = 

5 

®c 

= compressor efficiency = 

0.85 

eip 

= turbine efficiency = 

0.9 

Q* 

® energy parameter = 

2 


With these parameters, the basic isobaric combustion con- 
ventional gas turbine has an efficiency of 0. 256. The 
parameters of the augmented gas turbines are: 

Ab " combustor area 

Am - mixing area 

Ae " wave compression area 

There are two modes of firing of the detonation ducts which 
are discussed in Appendix A. In the single firing mode, 
the duct is operated at its nominal conditions, fuel is 
injected and the mixture is detonated. The flow in the 
duct is then allowed to return to its nominal conditions 
before it is detonated again. With 2 detonation ducts 
associated with each combustor, the maximum frequency is 
known to be controlled by the time for the detonation gases 
and associated wave systems to exit the duct (6* =3) and 
for flow to be re-established, which is also = 3 . Thus 
in Q =6, which is the shortest theoretical period of the 
system, the detonation ducts will fire twice. The other 
limit is obtained when the detonation ducts have no flow 
going through them. Calculated efficiencies for the single 
firing mode of operation are shown in Table 1. 

The detonation ducts can also be fired after the 
upstream contact surface exits the duct and the detonation 
duct volume is occupied by compressed gas. In principle, 
the duct could be fired the third time, but the pressure 
waves retard the flow significantly and the firing fre- 
quency diminishes with each firing quite rapidly. With 
double firing and two detonation ducts for each combustor, 
the minimum period is ^ = 10.2. Calculated values of 
efficiencies are shown in Table 2 . 
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It is clear that augmented gas turbines offer the 
potential for significant improvements over the conven- 
tional gas turbine efficiency of 0.256. Efficiency is 
higher for single firing mode of operation because greater 
fractions of the duct flow are processed by the detonation 
and shock waves. Efficiency also increases with decreasing 
mixer area, but the results for the lowest values of Am are 
probably not realizable because stable ejector operation 
may not be achievable. For Ab = 1 the combined detonation 
ducts and combustor exit area is 3.73 and a mixing area of 
4 is probably realistic. The effect of the wave compres- 
sion area on efficiency is minor so that it does not merit 
much consideration at this stage of the analysis. For Ab = 
2 the combined areas are 5.4 and a mixing area of 6 is 
probably quite reasonable. Representative results for a 
single firing mode are therefore: 

Ab = 1 V = 0.A5 

=2 = 0.37 

Similarly, for a double firing mode 

Ab = l T7=0.31 

= 2 = 0.29 

Clearly the single firing mode is superior and significant 
improvements over the conventional gas turbine efficiency 
of 0.26 are indicated. Even if the losses in the mixing 
and compression are greater than are currently estimated, 
the potential for enhancing the performance of low pressure 
ratio gas turbines has been demonstrated. 
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Table 1. Efficiencies of Augmented Gas Turbines 

Single Firing 


Ab 

Am 

Ae 

Xb 


1 

3 

3 

.350 

.256 

1 

3 

4 

.350 

. 534 

1 

3 

5 

.350 

.537 

1 

3 

6 

.350 

.539 

1 

4 

4 

. 350 

.452 

1 

4 

5 

. 350 

.459 

1 

4 

6 

.350 

.461 

2 

5 

5 

. 537 

.418 

2 

5 

6 

.537 

.422 

2 

5 

6 

.537 

( 

.373 ■ 

2 

6 

7 

.537 

.376 


Note - All areas are normalized by A<j 
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Table 2 . Efficiencies of Augmented Gas Turbines 

Double Firing 


Ab 

Am 

Ae 

Xb 


1 

2.5 

3.0 

.285 

.406 

1 

2.5 

4.0 

.285 

.411 

1 

2.5 

5.0 

.285 

.412 

1 

3.0 

3.0 

.285 

.368 

1 

3 . 0 

4.0 

.285 

.375 

1 

4 . 0 

4.0 

.285 

.313 

1 

4 . 0 

4.0 

.285 

.318 

2 

3 . 0 

4.0 

.461 

. 385 

2 

3.0 

o 

• 

in 

.461 

.387 

2 

3.0 

6.0 

.461 

.388 

2 

4.0 

4.0 

.461 

.332 

2 

4.0 

5.0 

.461 

.336 

2 

4.0 

6.0 

.461 

.338 

2 

o 

• 

in 

5.0 

. 461 

.293 

2 

5.0 

6.0 

.461 

.296 


Note - All areas normalized by A^ 
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ENGINE PERFORMANCE 


The basic calculations presented in the previous 
section indicated the potential of the detonation wave 
augmentation in terms of thermodynamic cycle efficiency. 
Here, the discussion will be in conventional terms of 
performance of various types of gas turbines. In the 
following, the representative optimum augmented gas turbine 
will be a single firing mode, Ab =1, Xb = 0.35, 17 = 0.45 
system. The off-design condition will be a conventional 
gas turbine with R = 5 and 17 = 0.26. With the conditions 
stated previously and propane as fuel, the combustor air 
fuel ratio is about 30. At design conditions, the detona- 
tion duct air/fuel ratio is about 45. The overall fuel/air 
ratio is therefore: 


Xb 1 - Xb 

mf/mair — + — 0.0261 

(mair/inf)b (niair/inf)b 

(19) 


Performance of a Pure Jet Engine 

The case considered is shown schematically in Figure 

11 . 



1 - INLET 

2 - COMPRESSOR 

3 - COMBUSTOR 

4 - DETONATION CHAMBERS 


5 - MIXING CHAMBER 

6 - TURBINE 

7 - NOZZLE 


Figure 11. Schematic Representation of a 
Detonation Wave Augmented Jet Engine 
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In this case, the net work is simply the increase in the 
kinetic energy of the gases, or: 


r) = (1 + Af) VE2/2qf 


( 20 ) 


With thrust F equal to the time rate of change of momentum, 
the fuel flow rate per unit of force is: 

1/2 

Wf/F = 3 (21) 

277qf(l + Af) 


In conventional terms, the performance of the optimum 
(maximum augmentation) and off-design engines is: 

Optimum wf = 0.83 Ib/lb hr 

Off-Design = 1.23 Ib/lb hr 

Here, qf for propane was taken as 22,000 Btu/lb. 


Performance of a Turbofan 


In this case, the net work is used to power the fan, 
as shown in Figure 12 . When the net work is equated to the 
change in kinetic energy of the total flow through the 
engine and the thrust is taken as the rate of change of 
momentum then: 

Ve 1 + 

Wf/F = (22) 

27jqf 1 + fjiX 

with 

M “ bypass ratio 
X - velocity ratio = Veb/"^E 

With eB = .85, X = 1/2, Ve = 2,200 ft/s, qf = 22,000 
Btu/lb, and fi- = 8 , the performance parameters are: 

Optimum Operation Wf/F = 0.345 Ib/lb hr 

Off-Design 
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0.596 



Figure 12 . 


Schematic Representation of the Turbofan Engine 


Performance of a Bypass Jet Engine 

In this configuration, the flow from the detonation 
ducts is used directly for propulsion as shown in Figure 
13. In this case, the turbine must deliver sufficient 
power to compress all the flow while only a fraction of it 
is going through the turbine itself. 

Compressor work 

Me Cp Ta(R^ - l)/ec 

Turbine work 

Mt Cp Tpf(l - R"^)et 

Equality of the above results in 

Mc/Mt = ®c ®t (Tbf/Ta) 
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The fuel input is 

Wf = [Mt/Aft +(Mc~Mt)/Afd]g 

with the assumption that the detonation duct flow is 
expanded to ambient pressure, the specific fuel consumption 
is: 


Wf/F 


1/Aft ( M -1 + 1/Afd) 



(23) 


here 


M = Mc/Mt 

xi = mass function of total flow 

In conventional units, the performance is calculated to be 

Optimum design wf/F = 2.4 Ib/lb hr 

Off-Design = 3.36 

The reason for the poor performance of this idea is -that 
over 60% of the flow from the detonation ducts is at a low 
temperature and pressure. The high pressure pulses, which 
in this case do not compress the rest of the flow, give 
sharp thrust pulses but of very short duration. 



Figure 13. Bypass Jet Engine 
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Performance of a Shaft Output Engine 


In this case, the whole output of the gas turbine is 
taken as shaft power. The arrangement is shown in Figure 
14. The specific fuel consumption is given in terms of the 
efficiency by the relation: 

Wf/W = 1/eqf (24) 

In conventional units of specific fuel consumption, the 
performance of the shaft engine is calculated to be: 

Optimum engine Wf/w = 0.257 Ib/hp hr 

Off-Design = 0.445 



Figure 14 . Schematic Representation of a 
Shaft Output Engine 


With the exception of the bypass jet engine, all the 
configurations indicate significant improvements in per- 
formance. Even if the full potential cannot be realized in 
a working engine, it appears very likely that detonation 
wave augmented engines could represent a major advance in 
technology. 


27 




CONCLUSIONS 


The study employed highly simplified, but realistic 
models of the gas dynamics of transient transverse detona- 
tion waves in transonic duct flows to estimate the maximum 
frequency at which such ducts could be operated. Concur- 
rently, a preliminary experimental effort not only proved 
the correctness of the hypothesis of the possibility of 
generating such waves, but also demonstrated that the 
induction distance for the formation of detonation waves 
was only about 1-2" and indicated that the pressure ratios 
across the wave were approximately lo. The combined 
theoretical and experimental results suggest that the 
salient features of the complex gas dynamic phenomena 
involved are modelled correctly. Therefore, the per- 
formance of detonation ducts operating as gas generators is 
predicted well enough for estimates of system performance. 
Losses due to mixing and sudden expansion into the mixing 
duct are accounted for using standard techniques, but it is 
possible that_ significant errors could result from the 
usual assumptions of c^asi-one dimensional flows. since 
the mixing process involves turbulent diffusion among 
streams at different pressures, calculated results of 
temperatures and velocities must be viewed with some 
caution. The pulse compression model for the transient 
flows of dissimilar volumes of gas contains the necessary 
conditions on energy and continuity equations but some 
uncertainty remains with regard to the ultimate state of 
the flow. The above considerations lead to the conclusion 
that even with highly reasonable models for the flow 
phenomena involved, considerable uncertainty remains as to 
the actual turbine entry conditions. 

Calculations of performance of the gas turbines 
operations in various configurations demonstrate signifi- 
cant enhancement of performance. At pressure ratios of 5 
the augmented gas turbines showed efficiencies of about 45% 
versus 26% for conventional engines. Even if the losses in 
the mixing and wave compression processes were doubled, a 
50% increase in efficiency still appears to be realizable. 
In view of the high cost of additional compressor stages 
which would be required to achieve such an improvement in 
performance, the probability of using very inexpensive duct 
remains should merit further attention. 
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The theoretical work in this field has gone as far as 
is practical. Far more elaborate flow models could be 
invented but the same uncertainties would still remain. It 
is concluded that the concept has been shown theoretically 
to have a high probability of success and must now be 
validated in a simple engineering laboratory experiment. 
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Af 

Cp 

Cv 

c 

D 

e 

E 
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g 

H 

Lo 

m 

M 

P 

qf 

Q 

Q* 

R 

t 

T 

Tmax 
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specific heat at constant pressure 

specific heat at constant volume 

acoustic speed 

detonation wave speed 

efficiency 

energy 

thrust 

gravitational constant 
heat release parameter 
length of the detonated volume 
mass 

Mach number 
pressure 

heating value of fuel 
heat 

heat release parameter 
compressor pressure ratio 
time 

temperature 

maximum temperature in the engine 
specific volume 
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Ve 

Wf 

w 

X 

Xb 

y 

r 

V 

A 

<P 

P 

d 

Subscripts 

a 

air 

b 

B 

bf 

c 

d 

D 

e 

f 

m 

o 


exit velocity 
fuel flow rate 
thermodynamic work 
axial distance 

fraction of flow through the combustor 
isentropic exponent 
1 - 1/y 

engine efficiency 
(3 - l)/2( y - 1) 

Tmax/Ta 

density 

dimensionless time 


ambient 

air 

burner or combustor 
bypass 

combustor exit conditions 

compressor 

detonated gas 

detonation wave 

wave compressor exit 

fuel 

mixing region 

detonated volume initial conditions 
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APPENDIX A 


GAS DYNAMICS OF THE 
DETONATION DUCT 


GAS DYNAMICS OF DETONATION DUCTS 


Presented here is an outline of the modelling and gas 
dynamics calculations which were performed as part of the 
feasibility study of detonation wave compression augmenta- 
tion of gas turbines. Surveys of open literature showed no 
comparable theoretical or experimental work so that this 
represents a new approach in the field of detonation waves. 
An independent laboratory experiment was performed to prove 
the existence of transverse detonation waves and to esti- 
mate some of the critical parameters. This work is des- 
cribed in Appendix B. Existing published work on detona- 
tion waves may be divided into three major categories: 

a. Studies of propagating steady waves 

b. Studies of stabilized waves in very high speed 
flows 

c. Studies of transitions from deflagrations to 
detonations 

The basic theory of steady detonation waves of the Chapman- 
Jouguet type has been known for about a century (e.g., 
Stanyukovichl) but numerous questions still remain. A 
detonation wave represents an essentially constant volume 
heat release with an associated increase in temperature and 
pressure. Extensive experimental studies of various 
explosive mixtures indicate that the fundamental Chapman- 
Jouguet hypothesis of sonic outflow from the detonation 
wave front is essentially correct. Calculations (e.g., 
Eisen, et al2) are generally in good agreement with experi- 
mental data. However, the nature of the detonation wave is 
itself in some question. While the Chapman- Jouguet detona- 
tion with sonic relative velocity of the exploded gas 
occurs most frequently, non-Chapman- Jouguet waves have also 
been observed (Gross and Oppenheim^, Schott^) . The latter 
cites instances of detonation products Mach numbers of 
1.1. Pressure ratios across detonation waves are also 
somewhat uncertain, with Eisen, et al2, calculating 

stoichiometric H 2 “Air pressure ratios to be 15.6, while 
JostS shows values of 15.2 and Sokolik^ shows a value of 
13. Similar uncertainties exist in the values for pres- 
sures and temperatures after the detonation wave is 
reflected from the wall as a shock wave. 
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When a detonation wave is stabilized in high speed 
flow, total pressure is diminished, but this is deemed 
acceptable in ramjets because intake losses are reduced. 
Such a loss could not be tolerated in gas turbines. 
Because of the continuing interest in high speed ramjets, 
the problem of stabilizing detonation waves in undiffused 
hypersonic flow has received considerable attention. Among 
the published studies are those of Waltrup, et al’7~9. The 
latter is an extensive survey of the field which updates 
that of DuggerlO. Numerous other experimental studies of 
standing detonation waves demonstrated the feasibility of 
stabilization of oblique detonation waves (e.g. , Nicholls 
and Daborall and Dunlap, et all2) and advantages of super- 
sonic combustion for certain ramjet applications were 
demonstrated by Dugger^^. while the field has been 
thoroughly explored, it is not clear that engineering 
applications for standing detonation waves do, in fact, 
exist. 

The problem of transition from deflagration to detona- 
tion in ducts filled with quiescent explosive mixtures has 
also received considerable attention. These studies have 
concentrated on the determination of the effects of size, 
surface roughness, residual turbulence, ignition mechanism, 
and the thermodynamic state of the explosive mixture. It 
is well known that wall roughness or obstacles which 
generate turbulence in the flow can reduce the transition 
distance by an order of magnitude (e.g., Sokolik^ and 
Brinkley and Lewisi^) . The configuration and the energy 
output of the igniter are also known to play an important 
role in the initial acceleration of the flame (Laderman, et 
all5) . Reference 14 describes experiments in which place- 
ment of obstructions near the igniter reduced the transi- 
tion or induction distance from 60 to 5 tube diameters, for 
a wide range of tube diameters. Reference 2 gives a 
transition distance for a quiescent ethylene-oxygen mixture 
of only 6 cm. Jost^ shows that detonations tend to form in 
about 3-4 tube diameters for tube diameters greater than 
about 15 mm. Jost also shows that the transition distance 
decreases by 15-40% when the initial temperature is 
increased form 15°C to 180°C. The speed of the detonation 
wave appears to be very weakly dependent on the mixture 
pressure and temperature. When the detonation is formed in 
short tubes or vessels then the precursor pressure waves 
can compress the unexploded mixture to a pressure ratio as 
high as 5 (Brinkley and Lewis^^) . a comprehensive 
parametric experimental study of Lee, et ali^ showed the 
effects of mixture composition and obstacles in the tube on 
induction distances in hydrogen-air mixtures. 
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All of the above studies are for detonation waves in 
stationary fluids and thus the results have very little 
direct bearing on the transverse waves in high speed flows 
which are of interest here. A comprehensive literature 
search using the extensive facilities of the UCLA Engineer- 
ing Library revealed no published data on transient trans- 
verse detonations in ducts. The only studies of the 
applications of travelling detonation waves in gas turbines 
appear to be those of Edwards!"^ who demonstrated that 
stable Chapman-Jouguet waves could be produced in an 
annulus. These results are interesting but they do not 
contribute to the gas dynamics problem considered here. 

The large number of uncertainties in the fundamental 
processes and details of the gas dynamics of the flow 
situation renders exact calculations essentially meaning- 
less and it is necessary to model the physics of the 
problem using reasonable approximations to the actual 
phenomena. Physical reasoning based on the experimental 
work performed in support of this program has led to the 
following mathematical models of the gas dynamics of 
transient, transverse detonation wave ducts. 


Initial Detonation 

Standing detonation waves were rejected for • gas 
turbine applications because of the very high total pres- 
sure losses associated with the Mach 4-6 shock waves. 
Initially, transient waves propagating longitudinally in 
ducts were given some consideration, but the idea was 
rejected because of the difficulties of uniform ignition 
across the duct which would be necessary for the establish- 
ment of a planar wave. Evaluation of numerous other ideas 
led to the conclusion that transient transverse waves were 
most likely to be successful. This was confirmed by the 
experiments described in Appendix B. The basic idea is to 
ignite a volume of detonable mixture from one or more sides 
of a rectangular duct or from the centerline of a circular 
duct. A line of symmetry or a solid surface are necessary 
to support the pressure which accelerates the combustion 
front and leads to the transition from deflagration to 
detonation. Schematic representations of the ignition 
models are shown in Figure A.l The model assumes that at 
the entrance to the detonation duct is a grid of fuel 
injectors which deliver fuel intermittently. Synchronized 
with the fuel injectors is the firing of the igniters at 
the sides or on the centerline of the duct. The latter is 


A-4 


not considered in depth at this time because of possible 
problems with unacceptably large induction distances for 
the formation of the detonation wave. After ignition the 
flame accelerates away from the wall while it is being swept 
downstream by the approximately Mach 1 channel flow. At 
some point a detonation wave is formed and is reflected 
from the far wall as a shock wave. The gas ahead of the 
detonation wave will be compressed by the precursor pres- 
sure waves generated by the accelerating flame. Some 
spreading of the flame in the longitudinal direction is 
inevitable as is turbulent diffusion so that the flame line 
will be quite irregular. It is also known (e.g., Jost^) 
that detonation waves tend to exhibit transverse instabili- 
ties which in tubes manifest themselves as spin. It almost 
certain that in the highly turbulent transonic duct flow 
some transverse instabilities will appear and a simple 
trajectory of the detonation wave is unlikely. The sudden 
increase of pressure in a volume of gas will generate shock 
waves in the longitudinal direction and the gas dynamic 
situation will be somewhat like that shown in Figure A. 2. 

It is known from the limited test data obtained from a 
simple detonation duct that the detonation wave forms 
within a short distance from the igniters (approximately 
1-2") . With duct width/length ratio in the ignition region 
of about 3-4 and a detonation wave Mach number of 4-6, the 
traverse time in the transverse direction is an order of 
magnitude lower than the longitudinal traverse time for an 
acoustic pulse. A reasonable engineering approximation to 
this very complex situation is therefore a concept of a 
sudden appearance of a volume of high pxessure and tempera- 
ture gas. This leads to a sudden expansion which generates 
shock waves which propagate away from the detonated gas. 
In the limit this model is that of a shock tube in which 
diaphragms are burst on both sides of the chamber of high 
pressure gas. 
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Figure A.l. Formation and Propagation of Detonation 
Waves Moving Transversely to the Flow 
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Figure A. 2. Expansion of the Detonation Products 
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The Shock Tube Problem 


When regions 1 and 2, initially at different tempera- 
tures and pressures, begin to interact with each other, a 
shock wave is propagated into the low pressure region and a 
train of rarefaction waves moves into the high pressure 
region. The two gases are separated by an interface at 
which pressure and velocities are equal. The situation is 
depicted in Figure A. 3 which also shows a shock wave 
entering an expanding duct. 



Figure A. 3. The Shock Tube Problem and Shock 
Wave in a Variable Area Duct 
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Analyses of the formation of the initial shock wave and the 
strength of the rarefaction waves follows the very clear 
exposition of Rudinger^S . The velocity of the gas in 
region 3 is given by the usual normal shock wave relation: 


Ue 


2aiMs(P3/Pl - 1) 

(A.l) 

y + 1 

( 7l - 1) (1 + .J: P 3 /P 1 ) 

7l - 1 


with a being the acoustic velocity, Ms the shock wave Mach 
number, P pressure, U velocity, 7 the isentropic exponent 
and subscripts denoting the corresponding regions. The 
Riemann invariant for region 4 is expressed as: 


U 4 = 


2 a 2 

^ 2-1 


1 - (P4/P2) 


y2zl' 

2 >'2 


(A. 2) 


with P 3 = P 4 and U 3 = U 4 


P3/P1 = P2/P1 


1 - 


^2 - 1 ai (P3/P1 - l)Ms 


2 ^2 
V2 - 1 


Tl - 1 a2 + 1 P3 

1 + 

'>'1 + 1 PlJ 


(A. 3) 


The relation between pressure ratio and shock Mach number 
is shown in Reference 19 to be: 


P 

Mg2 = ^ ^ + 1 _3 

2 7 2 7 Pj_ 


(A. 4) 


Solution of Equations A. 3 and A. 4 leads to the calculations 
of conditions across the shock and in the rarefaction fan. 
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Thus: 


T 3 /T 1 = (2/yi+l)2(yi Ms2 -(y-l]/ 2 )(l + (yi-l)/ 2 .Ms 2 )/Ms 2 


and 


(A. 5) 


P3/P1 


h + 1 

2 


Ms2 / (1 + 


y- 1 


2 



(A. 6 ) 


Conditions in region 4 are evaluated by means of isentropic 
relations. These solutions constitute the initial condi- 
tions for the interaction phenomena which are studied using 
the method of characteristics. 

The representative calculations exhibited in the 
following are based on the following initial conditions: 


yi =1.4 

Y2 = 1.25 

^2/^1 - 2.5 

P2/P1 = 8.5 


These values are based on theoretical considerations- and 
preliminary test data. The initial values corresponding to 
these parameters are: 


Ms = 

1.99 

as/ao 

■ = 1.244 


P3/P0 = 

4.47 

V^o = 2.66 

T3/T0 = 

1.68 

In all of the 
in every flow 

following 

region 

calculations, it was assumed 
the gas was perfect and that 

that 

the 


isentropic exponents were different but constant. 
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Initial Interactions 


The initial expansion of the detonated mixture creates 
rarefaction wave fans, which intersect within the detonated 
volume, as shown in Figure A. 4. The numbers on the char- 
acteristics are the Riemann invariants which in this case 
are defined as: 


I = 500[1 + a/ao + ~ ^ U/ao] (A. 7a) 

2 

II = 500[1 + a/ao - ^ U/ao] (A. 7b) 

2 

This identification of the invariants is purely arbitrary 
and is used here simply for convenience. The fundamental 
statement is that the quantities [a + ( 7' - l)U/2] remain 

constant on trajectories defined by U + a. 

The changes in state variables indicated in the above 
Figures are due to the interactions of the initial and 
reflected rarefaction wave trains and the passage of 
interfaces and shock waves. 

The dimensionless time for the problem is: 


d — tao/Lo 


and the dimensionless coordinate is X/Lo with X measured 
from the downstream boundary of the detonated volume. 
Variation of pressure with time throughout the detonated 
volume is shown in Figure A. 5. The same variations at 
later times are shown in Figure A. 6. All of the above 
quantities are shown in a coordinate system which moves at 
Mach 1 with the flow in the duct. 
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DIMENSIONLESS TIME - 9 = ta^/L 



Figure A. 4. Characteristics Net for the Initial 
Interaction in a Moving Coordinate System 
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Figure A. 5. Distribution of Pressure in the Detonation 
Volume in a Moving Coordinate System 
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PRESSURE RATIO - P/P 



AXIAL DISTANCE - X/L^ 


Figure A. 6. Distributions of Temperature and Pressure 

within the Detonation Gas 
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TEMPERATURE RATIO - T/T 




Motion of Discontinuities 


After the rarefaction wave trains intersect in the 
detonated volume, the waves intercept the contact surfaces. 
Interception -of a contact surface by a rarefaction wave 
results in the deceleration of the contact surface, reflec- 
tion of a rarefaction wave and a transmission of a weak 
rarefaction wave. The situation is shown in Figure A. 7. 


t 


X 



// 

// 






/ 


Figure A. 7. Interaction of a Rarefaction Wave 
with a Contact Surface 


Conditions for the reflected and transmitted waves are 
obtained from considerations of the Riemann invariants and 
the continuity of pressure and velocity at the contact 
surface. The calculated results for the initial interac- 
tion of the rarefaction waves with the contact surface are 
shown in Figure A. 8. The transmitted waves propagate at 
sonic velocity and therefore overtake the shock wave. The 
resulting flow situation is somewhat similar to that shown 
in Figure A. 7 except that there is no transmitted wave. 
Calculations of interactions of waves and contact surfaces 
in the inertial coordinate system of the duct are shown in 
Figure A. 9. During the interactions, the flow is being 
swept downstream (to the left) at Mach 1 so that the up- 
stream extent of the gas dynamic phenomena is rather 
limited. It should be noted that the contact surface which 
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initially moved upstream is swept back to its initial 
location at 0 = 1.1 and clears the detonated volume at 
$ - 2 . 1 . Therefore, in principle, the compressed gas in 

the detonation volume could be mixed with a new charge of 
fuel and detonated at some time 6 :^ 2 . 1 . This is one of 
the alternative modes of operation which is considered 
later. 

The shock wave which has an initial Mach number of 
about 2 , in fact initially propagates upstream in the 
inertial coordinate system at a speed corresponding to Mach 
1 since the duct flow is taken to be to the left at Mach 1. 
The first rarefaction waves overtake the shock wave at the 
point X/Lo - 2.1 or about 1 detonated volume length 

upstream of the igniters. Motions of the shock wave and 
the upstream contact surface are shown in Figure A. 10. It 
should be noted that at a distance of 2.4 L© from the 
igniters, the shock wave Mach number is reduced to about 
1.4. Correspondingly, the pressure ratio has diminished 
from 4.5 to 2.1. This indicates that the rarefaction waves 
from the opposite side of the detonated volume alternate 
the upstream propagating shock wave quite rapidly. The 
phenomena of shock wave and contact surface interactions 
downstream of the detonated volume will be considered in 
the discussion of the duct output. 
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DIMENSIONLESS TIME - 6> = ta /L 
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Figure A. 8. Interactions with the Detonated Gas and the 
Shock Wave in a Moving Coordinate System 
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DIMENSIONLESS DISTANCE - X/L^ 

Figure A. 9. Salient Features of the Primary Detonation 

in Inertial Coordinates 
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Attenuation of Upstream Shock Waves 


The shock waves propagating upstream from the deto- 
nated volume must be attenuated to minor acoustic distur- 
bances in order to minimize the pressure fluctuations which 
might affect compressor performance. It is known that 
rapid area increases are very effective in reducing the 
strength of transmitted shock waves. In general, the gas 
dynamic situation is as shown in Figure A. 11. 



Figure A. 11. Interactions of Shock Waves 
with Area Changes 


Here, following Rudinger^®, the area change is shown to be 
instantaneous. This approximation is discussed by Refer- 
ence 18 and is shown to give good agreement with experi- 
mental results. The more elaborate analysis of Whitham^O 
which considers gradual area changes could be adapted to 
the solution of shock waves moving through variable area 
ducts with gas flows, but since only an estimate is needed 
here, the considerable analytical effort involved could not 
be justified at this point. 

The actual flow situation will vary significantly with 
the strength of the incident shock wave and the Mach num- 
bers of the flow in sections 1 and 2 . Solutions of the 
problem are obtained by matching the conditions in region 5 
which results from an isentropic expansion through the 
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reflected waves with isentropic flow between area changes 
from region 5 to region 6 and in turn matching the velocity 
and pressure at the contact surface between regions 6 and 
4 . The process is quite straightforward and the required 
large number of interactions is easily performed on a com- 
puter. The procedure is as follows: 

a. The velocity of the gas in region 3 is computed 
using: 

U3/ai = Ml - Ms(l - Ri3) (A. 8) 

R13 = ^1/^3 

b. The temperature ratio T 3 /T 1 is calculated using 
shock wave relations 

c. Mach number in region 3 is calculated using 

U 

M 3 = _ 2 (Ti/T 3 ) 1/2 (A. 9) 

a 

1 

d. With a chosen value of P 5 the Mach numbers in 
region 3 and 5 are connected using isentropic flow 
relations 

e. Velocity in region 5 is calculated using isentro- 
pic flow relations 

f. Transmitted shock wave Mach number is guessed and 
conditions in region 4 are calculated 

g. Using isentropic relations and conditions in 
region 5, the conditions in region 6 are calcu- 
lated 

h. Calculated conditions in regions 6 and 4 are com- 
pared and the guesses for Mt and P 5 are adjusted. 


The process is repeated until the contact surface 
conditions are satisfied to an arbitrary degree of preci- 
sion. Solutions for strong incident shock waves in air 
flowing at Mi = 1 are shown in Figure A. 12. 


A-20 



AREA RATIO 
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Figure A. 12. Mach numbers of the Shock Waves Transmitted 

Out of the Detonation Duct - 


In the particular situation considered here, the 
incident shock wave Mach number is 1.4 so that the problem 
is simplified greatly. Since the shock wave is moving into 
a counterflow at M = 1, its net speed in inertial coordi- 
nate system is only M = 0.4 and the gases in region 3 are 
swept downstream. The contact surface also moves down- 
stream and the transmitted shock wave is extremely weak. 
In fact, at an area ratio of 1.34 which corresponds to a 
flow Mach number of 0.5, Mt =1.02 and the pressure ratio 
is only 1.047. At an area ratio of 2, which corresponds to 
a flow Mach number of 0.3, the transmitted shock wave 
decays to an acoustic pulse. 

I 

Detailed calculations employing inviscid flow assump- 
tions are not very meaningful when the shock waves are 
weak because of strong dissipation effects. However, the 
results indicate that even if different fuels are used, or 
the distance between the duct inlet and the detonated 
volume is reduced, shock waves with Mach numbers of 1.5 
2 , 0 can be attenuated by area ratios of less than 5 . 
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Output of the Detonation Duct 


The detonation duct is a gas generator whose outpur 
varies with time. Here, and in all results exhibited in 
this Appendix, it is assumed that the air in the duct flows 
at sonic velocity, the compressor pressure ratio is 5:1, 
and that a stoichiometric mixture of propane and air is 
detonated. It is shown in Figure A. 9 that the upstream 
contact surface clears the detonated volume after i9 = 2.1 
and the flow is fairly uniform at (9 = 3. Similar situation 
exists at 9 = 8. In all of the exhibited results, it will 
be assumed that the duct is fired at 6 =0, 3, 8. After 
each firing, the flow velocity in the duct is diminished 
and the cyclic frequency is reduced. 

Summary of the thermodynamic states and flow condi- 
tions at the duct exit, which is defined as X/Lq = -1 is 
presented in Table A.l for 3 consecutive firings. Obvi- 
ously, in actual operation the cut could be fired 1, 2, or 
3 times. Variations of the static pressures, temperatures, 
and Mach numbers with dimensionless time are shown in 
Figure A. 13. All quantities are normalized by the initial 
duct flow parameters which are denoted by the subscript o. 
Corresponding total pressures and temperatures are shown in 
Figure A. 14. Flow velocities at the duct exit are shown 
in Figure A. 15. 

The most important parameter for the gas generator 
aspect of the detonation duct is the power output. Since 
only comparable quantities are of interest at this stage, 
the instantaneous power output is defined as isentropic 
expansion through a turbine to the initial pressure at the 
compressor inlet. This output is normalized by the energy 
flow rate through the detonation duct 

• 

Wq = Pq Uq Cp To (A. 10) 


The output from such isentropic blowdown turbines supplied 
by the detonation ducts is compared with ideal expansion 
from a combustor with flow at Mq = 0.25 in Figure A. 16. 
Comparisons of the total energy generated in a complete 
cycle consisting of 3 firings and recovery of the flow to 
initial conditions is 0.4 shown in Figure A. 17. Obviously, 
the highly simplified calculations shown here must be 
viewed as indications of the potential of the detonation 
duct as a gas generator rather than from engineering quan- 
tities . 
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Table A.l. Summary of Duct Output 


EVENT 


P/Po. 

T/To 

pt/Po 

Tt/To 

U/ao 

M 




FIRST 

FIRING 




0 

0.00 

1.20 

1.00 

1.89 

1.20 

1.00 

1.00 

A 

0.33 

4.47 

1.64 

23.90 

2.65 

2.24 

1.75 

B 

0.46 

4.47 

5.48 

7.26 

6.06 

2.27 

0.97 

C 

0.49 

4.47 

5.48 

7.26 

6.13 

2 . 24 

0.96 

D 

0.51 

4.20 

5.43 

6.88 

5.99 

2.12 

0.91 

E 

0.53 

3.90 

5.35 

6.08 

5.84 

2 . 00 

0.86 

F 

0.59 

3.20 

5.14 

4.36 

5.47 

1.62 

0.72 

G 

0.73 

2.26 

4.80 

2.57 

4 . 92 

1.00 

0.46 

H 

1.25 

2.20 

4.77 

2.49 

4 . 89 

0.96 

0.44 

I 

1.72 

2 . 11 

4.73 

2.34 

4.83 

0.89 

0.41 

J 

2.82 

2.11 

4.73 

2.34 

4.83 

0.89 

0.41 

K 

2.88 

1.22 

1.12 

1.97 

1.27 

0.86 

0.81 




SECOND 

1 FIRING 




0 

0.00 

1.22 

1.12 

1.97 

1.27 

0.86 

0.81 

A 

0.35 

5.46 

1.81 

22.10 

2.45 

2.11 

1.57 

B 

0.50 

5.46 

6.06 

8.47 

6.62 

2.11 

0.86 

C 

0.50 

5.46 

6.06 

8.47 

6.62 

2.11 

0.86 

D 

0.51 

5.13 

5.99 

7.59 

6.47 

1.98 ■ 

• 0.8-1 

E 

0.53 

4.76 

5.90 

6.74 

6.32 

1.85 

0.76 

F 

0.60 

3.91 

5.67 

4.91 

5.93 

1.45 

0 . 61 

G 

0.75 

2.76 

5.29 

2.97 

5.37 

0.80 

0.35 

H 

1.27 

2.70 

5.26 

2.81 

5-. 34 

0.76 

0.33 

I 

1.77 

2.57 

5.22 

2.55 

5.25 

0.68 

0.30 

J 

4.21 

2.57 

5.22 

2.55 

5.25 

0.68 

0.30 

K 

4.22 

1.49 

1.23 

1.91 

1.32 

0.67 

0.61 




THIRD 

FIRING 




0 

0.00 

1.49 

1.23 

1.57 

1.32 

0.67 

0.61 

A 

0.32 

6.66 

2.00 

20.85 

2.77 

1.96 

1.39 

B 

0.46 

6.66 

6.70 

9.13 

7 . 15 

1.87 

0.72 

C 

0.46 

6.66 

6.70 

9.13 

7.15 

1.87 

0.72 

D 

0.47 

6.25 

6.62 

8.49 

7.04 

1.87 

0.71 

E 

0.48 

5.80 

6.52 

7.57 

6.88 

1.69 

0.66 

F 

0.56 

4.77 

6.27 

5.59 

6.48 

1.28 

0.51 

G 

0.70 

3.37 

5.85 

3.91 

6.03 

0.59 

0.24 

H 

1.20 

3.29 

5.82 

3.40 

5.86 

0.56 

0.23 

I 

1.66 

3.14 

5.77 

3.22 

5.80 

0.46 

0.19 

J 

5.66 

3 . 14 

5.77 

3.22 

5.80 

0.46 

0.19 

K 

5.80 

1.82 

1.38 

2.13 

1.44 

0.56 

0.48 
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KEY TO TABLE A.l 


0 Initiation of detonation 

A Shock crossing 

B Contact surface crossing 

C First rarefaction wave 

D Second rarefaction wave 

E Third rarefaction wave 

F Middle rarefaction wave 

G Final rarefaction wave 

H First wave reflected from the contact surface 

1 Final wave reflected from the contact surface 

J Contact surface crossing 

K Reflected wave from the shock 

More realistic power extraction models are used in the 
main report to evaluate potential gas turbine performance 
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DIMENSIONLESS TIME - 6 = ta„/L 



Figure A. 13. Variations with Time of Mach Number, 
Static Temperature and Static Pressure 
at the Duct Exit 
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DIMENSIONI ESS TIME - 6= tajl 



0 5 10 15 20 25 


TOTAL TEMPERATURE AND PRESSURE 


Figure A. 14. Gas Dynamic Conditions at 
Duct Exit (X/Lq = “1) 
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DIMENSIONLESS TIME - $ = ta„/L 



Figure A. 15. Variation with Time of Gas Velocity 

at Duct Exit 
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DIMENSIONLESS TIME 




Figurs A. 17. Energy Generated by the Detonation Duct 
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f I 




List of Symbols 


a Acoustic velocity 

Cp Specific heat at constant pressure 

I Characteristic value 

II Characteristic value 

Lq Length of the detonated volume 
M Mach number 

P Pressure 

Rij Density ratio = ^i/Pj 

t Time 

T Temperature 

U Particle velocity 

W Thermodynamic work or energy 

X Axial distance 

y Isentropic exponent of the gas 

p gas density 

0 Dimensionless time 


Subscripts 

o 

1 , 2 , • • . 
c 
s 
T 


Initial conditions before detonation 

Regions 

Combustor 

Shock 

Transmitted shock 
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APPENDIX B 


EXPERIMENTAL STUDIES OF 
DETONATION DUCTS 


EXPERIMENTAL STUDIES OF DETONATION DUCTS 


A simple experimental program was carried out in 
support of, but not as part of, the contract. Due to 
limitations of time and finances, the experimental program 
was aimed at demonstrating that transverse detonation waves 
could be produced, and to indicate the pressure levels 
which could be achieved. 

From the beginning of the program of studies of 
detonation wave augmented gas turbines, serious questions 
were raised regarding the possibility of forming transverse 
detonation waves in transonic flows. With the exception of 
the rotating detonation waves of Edwards, which was noted 
in Appendix A, all previous work in detonation waves had 
been done in quiescent gases or stabilized waves in hyper- 
sonic flow. Because of the range of occasionally conflict- 
ing experimental data on the initiation and propagation of 
detonation waves, it became obvious that some experimental 
data were needed to support the fundamental assumptions 
underlying the theoretical work. 

A rectangular detonation duct was fabricated out of 
heavy steel plates. Details of the duct and the assembled 
equipment are shown in Figures B.1-B.4. . Pressures -were 
measured using Kistler 602A high frequency pressure trans- 
ducers mounted in instrumentation plugs such as may be seen 
in Figure B.5. Air was supplied continuously from a high 
pressure wind tunnel reservoir. Propane was injected 
intermittently from various tanks. Ignition was found to 
be the easiest and most effective when a row of spark plugs 
was fired by a series of high frequency interrupters or an 
automotive distributor. There was no direct control over 
the air/fuel ratio so that various air supply and propane 
reservoir pressures had to be tried until detonations could 
be obtained on every attempt. Occurrence of detonations 
was determined from transducer pressure outputs, and also 
audibly from the very loud sharp cracks, rather than dull 
thuds which followed a misfire or delayed ignition. 

Representative results are shown in Table B.l. The 
relatively reliable velocities which are derived from the 
time increments between the pressure recording at two 
transducers correspond to shock waves whose pressure ratios 
are much higher than those exhibited in the Table. There 
was some uncertainty regarding the pressure level in the 
duct and it is possible that the detected pressures at 
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different transducers may have been caused by different 
waves. The data does indicate that detonation waves are 
formed in very short distances from the igniters (about 
1 - 2 ") . 

In the interest of using conservative basic data in 
the theoretical work, the bulk of the pressure data, which 
indicated pressure ratios of 9.6, was taken as the basis 
and reduced to 8.5 in the calculation. 
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13.500 6 EQUAL SPACES AT 

3.500 * 21.000 
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4 MATERIAL : MILD STEEL . 

3 DIMENSIONAL TOLERANCES: t.XXX*t.OIO, i.XX* + .03 

2. DIMENSIONS ARE IN INCHES . 

I. REMOVE ALL BURRS AND SHARP CORNERS . 

NOTES : UNLESS OTHERWISE SPECIFIED 



■BATTERY 


RELAY 


12.000 


© 


3.500 MAX 


■SPARK PLUG WIRE 


-TRANSDUCER 


51.000 


6 

TRANSDUCER 
RETAINING BOLT 

5 

HEX HEAD BOLT 

4 

LEFT SIDE 
PLATE 

3 

RIGHT SIDE 
PLATE 

2 

BOTTOM PLATE 

1 

TOP PLATE 

ITEM 

NO. 

DESCRIPTION 


STEEL 


MILD STEEL 


MATERIAL 


Figure 4.B. Detonation Duct Assembly 
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Figure 5.B. Details of an Instrumentation Plug 
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